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Abstract 

The scaling form of the critical equation of state is computed for 0(iV)-symmetric models. 
We employ a method based on an exact flow equation for a coarse grained free energy. A 
suitable truncation is solved numerically. 
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A precise computation of the critical equation of state near a second order phase transition 
is an old problem. From a general renormalisation group analysis 0 one can prove the Widom 
scaling form FT = cj)^f{{T — Tc)/ct>^^^) for the relation between the magnetic field H, the 
magnetisation (f> and the difference from the critical temperature T — Tc- In several models the 
critical exponents /3 and 5 have been computed with high accuracy Q but the scaling function / 
is more difficult to access. Previous attempts include an expansion in 4 — e dimensions in second 
order in e (third order for the Ising model) |^. A particular difficulty for a direct computation 
in three dimensions arises from the existence of massless Goldstone modes in the phase with 
spontaneous symmetry breaking for models with continuous symmetry (e.g. Heisenberg models 
with 0{N) symmetry for > 1). They introduce severe infrared problems within perturbative 
or loop expansions. 

Recently a non-perturbative method has been proposed which can systematically deal with 
infrared problems. It is based on the average action F^ [Q| which is a coarse grained free energy 
with an infrared cutoff. More precisely F^ includes the effects of all fluctuations with momenta 
but not those with q'^ < k‘^. In the limit /c —> 0 the average action becomes the standard 
effective action (the generating functional of the IPI Green functions), while for A: —> oo it equals 
the classical or microscopic action. It is formulated in continuous space and all symmetries of 
the model are preserved. There is a simple functional integral representation Q of F^ also for 
k > 0 such that its couplings can, in principle, also be estimated by alternative methods. 

The exact non-perturbative flow equation Q for the scale dependence of F^ takes the simple 
form of a renormalisation group improved one-loop equation |^] 
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The trace involves a momentum integration and summation over internal indices. Most impor¬ 
tantly, the relevant infrared properties appear directly in the form of the exact inverse average 
( 2 ) 

propagator F^ , which is the matrix of second functional derivatives with respect to the fields. 
There is always only one momentum integration - multi-loops are not needed - which is, for a 
suitable cutoff function Rk{q^) (with F?fc(0) ~ kP', Rki^f oo) ~ both infrared and 

ultraviolet finite. 

The flow equation i) is a functional differential equation and an approximate solution 
requires a truncation. Our truncation is the lowest order in a systematic derivative expansion 
ofFfc iga 






( 2 ) 


Here (j)^ denotes the iV-component real scalar field and p = We keep for the potential 

term the most general 0(A^)-symmetric form Uk{p) since Uo{p) encodes the equation of state. 
The wavefunction renormalisation is approximated by one A:-dependent parameter Z^- Next 
order in the derivative expansion would be the generalization to a p-dependent wavefunction 
renormalisation Zk{p) plus a function Yk{p) accounting for a possible different index structure 
of the kinetic term for Ai > 2 [Q, Going further would require the consideration of terms 
with four derivatives and so on. Goncerning the equation of state for the present model, the 
omission of higher derivative terms in the average action typically generates an uncertainty of 
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the order of the anomalous dimension rj. The main reason is that for 77 = 0 the kinetic term 
in the fc-dependent inverse propagator must be exactly proportional to both for —> 0 and 

^ 00 . For the three-dimensional scalar theory r/ is known to be small and the derivative 
expansion is, therefore, expected to give a reliable approximation. This holds for arbitrary 
constant “background” field (/>“. Similar, although less stringent, arguments indicate a weak 
p-dependence of the kinetic term. For the scaling solution for N = 1 this weak p-dependence 
has been established explicitly Q. 

In this letter we compute the effective potential (Helmholtz free energy) lim^^o Uk{p) = U (p) 
for the 0(A^)-model directly in three dimensions from a solution of eqs. dU), (H). We extract the 
Widom scaling form of the equation of state and give semi-analytical expressions for N = 1 and 
N = 3. Its asymptotic behavior yields the universal critical exponents and amplitude ratios. An 
alternative parametrisation of the equation of state in terms of renormalised quantities is used 
in order to compute universal couplings. 

For a study of the behavior in the vicinity of the phase transition it is convenient to work 
with dimensionless renormalised fields 


p = Zkk V, Uk{p) = k ^Uk{p{p)). 
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With the truncation of eq. (§) the exact evolution equation for = duk/dp ^ reduces to 
the partial differential equation 


^ = (-2 + p) 4 + (1 + p) pu'i 


~ ^ (3'Wfc + ‘^pu'k'j ll [u'k + 2pu'^; p), (4) 

where t = ln(A:/A), with A the ultraviolet cutoff of the theory. The anomalous dimension p is 
given in our truncation by 


P = 


d 2.0 

InZfc = ■^KA^m2,2(2A«;). 
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with K, the location of the minimum of the potential, = 0, and A the quartic coupling, 

u'1{k) = A. The “threshold” functions l\ and 2 result from the momentum integration on 
the r.h.s. of eq. (^) and account for the decoupling of modes with effective mass larger than k. 
They equal constants of order one for vanishing arguments and decay fast for arguments much 
larger than one. For the choice of the cutoff function R^, employed here their explicit form can 
be found in refs. |^, ^]. 

To obtain the equation of state one has to solve the partial differential equation (^) for 
A; —> 0. Algorithms adapted to the numerical solution of eq. have been developed previously 
1^] and we refer to this work for details. The integration starts at some short distance scale 
k~^ = A“^ {t = 0) where the average potential is equal to the microscopic or classical potential 
(no integration of fluctuations has been performed). We start with a quartic classical potential 
parametrized as u'f^{p) = X\{p—K\). In the phase with spontaneous symmetry breaking the order 
parameter po = lim^^o-^^^^^ takes a non-vanishing value. In the symmetric phase the order 
parameter vanishes, i.e. po = 0 for /c = 0. The two phases are separated by a scaling solution for 
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which u'^{p) becomes independent of k. For any given Aa there is a critical value Kcr for which 
the evolution leads to the scaling solution. A measure of the distance from the phase transition 
is the difference 5k,\ = k\ — Kcv If is interpreted as a function of temperature, the deviation 
5k\ is proportional to the deviation from the critical temperature, i.e. = A{T)(Tc — T) with 

A{n) > 0 . 

The external field H is related to the derivative of the effective potential U' = dU/dp by 
Ha = U'(f)a- The critical equation of state relating the temperature, the external field and the 
order parameter can then be written in the scaling form ((f) = -v/ 2 p) 


^5-1 HX), X 
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with critical exponents 5 and /3. For ^ > oo our numerical solution for U' obeys U' ~ (f)^~^ 

to high accuracy. The inferred value of 6 is displayed in the table, and we have checked the 
scaling relation 6 = (5 — rj)/(1 + p). The value of the critical exponent p is obtained from eq. (^) 
for the scaling solution p. We have also verified explicitly that / depends only on the scaling 
variable x for the value of (3 given in the table. In figs. ^ and |2| we plot log(/) and \og(df /dx) 
as a function of log|x| for N = 1 and N = 3. Fig. corresponds to the symmetric phase (x > 0) 
and fig. to the phase with spontaneous symmetry breaking (x < 0 ). 

One can easily extract the asymptotic behavior from the logarithmic plots. The curves 
become constant both for x —> O"*" and x —> 0“ with the same value, consistently with the 
regularity of /(x) at x = 0. For the universal function one obtains 


lim f{x) = D (7) 

x^O 

and H = on the critical isotherm. For x —> oo one observes that log(/) becomes a linear 
function of log(x) with constant slope 7 . In this limit the universal function takes the form 


lim /(x) = {C~^) ^x"^, 

X^OO 


( 8 ) 


or lim^^o = m^, and we have verified the scaling relation 7 //? = 

5 — 1 . One observes that the zero-field magnetic susceptibility, or equivalently the inverse 
unrenormalised squared mass m~‘^ = x, is non-analytic for —> 0 in the symmetric phase: 

X = In this phase the correlation length ^ = (Zqx)^'^^) which is equal to the inverse 

of the renormalised mass mu, behaves as ^ with v = p/(2 — p). 

In the phase with spontaneous symmetry breaking (x < 0) the plot of log(/) fig. shows a 
singularity for x = i.e. 

f(x = = 0. (9) 

The order parameter for H = 0 therefore behaves as (f> = B(5ntC)^. Below the critical tempera¬ 
ture the longitudinal and transversal magnetic susceptibilities xl and xt are different for N > 1 
(/' = df/dx) 


Xl = 


d‘^U 






( 10 ) 
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This is related to the existence of massless Goldstone modes in the (N — 1) transverse directions 
which imply that the transversal susceptibility diverges for vanishing external field. Fluctuations 
of these massless modes also induce a divergence of the zero-field longitudinal susceptibility. This 
can be seen from the singularity of the plot of log(/^) for iV = 3 in fig. The first derivative 
of the universal function with respect to x vanishes as > 0, i.e. f'{x = —= 0 
for N > 2. For = 1 there is a non-vanishing constant value for f'{x = with a 

finite zero-field susceptibility x = C'“(5 ka)~"’' where (3. For a 

non-vanishing physical infrared cutoff k the longitudinal susceptibility remains finite also for 
N >2: XL (kpo) In the ordered phase the correlation length for iV = 1 behaves as 
^ and, also for > 1, the renormalised minimum pqr = Zqpq of the potential U 

scales as poR = E{5 k\Y. 

The amplitudes of singularities near the phase transition D, B and E are shown in 

the table. They are not universal since different short distance physics will result in different 
wavefunction renormalisations and Z^p 2 . All models in the same universality class can, 
however, be related by a multiplicative rescaling of (j) and 5ka (or Tc — T) resulting in x ^ CxX 
and / ^ Off. Ratios of amplitudes which are invariant under this rescaling are universal. We 
display the universal combinations jC~ ^ ^ ~ DB^~^, 

and E in the table. 

The asymptotic behavior observed for the universal function can be used to obtain a semi- 
analytical expression for f{x). We find the following fit to reproduce the numerical values for 
both / and df /dx within 1 % deviation (apart from the immediate vicinity of the zero of / for 
N = 3, cf. eq. (|I3)): 

ffitix) = D(l + B^/'^x^ (l + 0 ®) , ( 11 ) 

with c = The parameter a is determined by the order of the pole of 

f~^ at X = —i.e. a = 1 (a = 2) for A^ = 1 (N > 1). The fitting parameters are chosen as 
0 = 0.569 (1.312) and A = 0.180 (-0.595) for A = 1 (3). 

There is an alternative parametrisation of the equation of state in terms of renormalised 
quantities. In the symmetric phase {6k\ < 0) we consider the dimensionless quantity 






m 
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s = —= -(e+)3(C+)-ix-2/3 
rriR 2 
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with PR = Zqp and = Zg The derivatives of F at s = 0 yield the universal couplings 


ds m-R niR ds^ 

and similarly for higher derivatives. They determine the behavior of / for x S> 1/2 


(13) 


fix) = (C+)-lx^ + i^(^+)3(C+)-2x^-2/3 + iz.«(^+)6(C+)-3x^-^^ + . . . (14) 

2 niR 8 


In the ordered phase {Skr > 0) we consider the ratio 

Gis) = ^ = ^f2f-3(-x)-^/(x), s = -^ = B-\-x)-^d_ 


POR 


POR 


(15) 
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The values for the universal couplings 


■^(1) = U'/Kpor) = z>ij 


(16) 


_ UniPoR) - 

ds poR poR ’ 


as well as XR/mR and ur are given in the table. One observes that for > 1 the renormalised 
quartic coupling Xr vanishes in the ordered phase. This results from the presence of massless 
fluctuations. For x near the scaling function is approximated by 


fix) = E^B-\-x)^(i-x)-‘^d _ ^ 2 ^ (2B^^ + VR(i-x)-^^ - B"^)) +... ( 17 ) 

In summary, our numerical solution of eq. gives a very detailed picture of the critical 
equation of state. The numerical uncertainties are estimated by comparison of results obtained 
through two independent integration algorithms [^. They are small, typically less than 0.3% 
for critical exponents and 1 — 3% for amplitudes. The scaling relations between the critical 
exponents are fulfilled within a deviation of 2 x 10“^. The dominant quantitative error stems 
from the truncation of the exact flow equation and is related to the size of the anomalous 
dimension rj ~ 4%. This is consistent with the fact that the critical exponents and amplitudes 
calculated here typically deviate by a few percent from the more precise values obtained by 
other methods p. If the equation of state is needed with a higher accuracy one has to extend 
the truncation beyond the level of the present work. 
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Tables 


N 

P 

7 

5 

u 

V 

^R/niR 


^r/por 

i'R 


1 

0.336 

1.258 

4.75 

0.643 

0.044 

9.69 

108 

61.6 

107 


3 

0.388 

1.465 

4.78 

0.747 

0.038 

7.45 

57.4 

0 

~ 250 



(T+ 

D 

B 


E 

C+/C- 




i+E 

1 

0.0742 

15.88 

1.087 

0.257 

0.652 

4.29 

1.86 

1.61 

0.865 

0.168 

3 

0.0743 

8.02 

1.180 

0.263 

0.746 

- 

- 

1.11 

0.845 

0.196 


Table 1: Parameters for the equation of state. 


Figures 

Fig. 1 ; Logarithmic plot of / and df /dx for x > 0. 
Fig. 2 ; Logarithmic plot of / and df /dx for x < 0. 







































































